X-ray scattering as a probe for warm dense mixtures and high-pressure miscibihty 
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We demonstrate the abilities ol elastic x-ray scattering to yield information on dense matter with 
multiple ion species and on the microscopic mixing in dense materials. Based on partial structure 
factors from ab initio simulations, a novel approach for the elastic scattering feature is applied to 
dense hydrogen-beryllium and hydrogen-helium mixtures. The scattering signal differs significantly 
between single species, real microscopic mixtures, and two separate fluids in the scattering volume. 
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Most materials in nature and technical applications 
consist of multiple atomic species and their properties 
can be only understood if all mutual interactions are con- 
sidered. This is also true for high energy density science 
applied to astrophysics and inertial confinement fusion. 
Indeed, objects such as old stars and giant gas planets 
can be considered as a natural laboratory for the equation 
of state of dense, high-pressure mixtures [l|-|Jl- Demixing 
and the subsequent segregation of the heavy elements can 
play a crucial role in the energy balance of these objects 
[5|. Moreover, the assumed separation into helium- rich 
and helium-poor phases strongly influences the evolution 
and present internal structure of gas planets and might be 
the reason for the observed low helium concentration in 
the atmosphere and the high luminosity of Jupiter and 
Saturn. Phase separation in hydrogen-helium mixtures 
has been predicted by quantum simulation for conditions 
found in the interior of gas giants [6, 7]. However, the 
verification of these theoretical predictions in an astro- 
physical setting is quite indirect and requires full scale 
planet modeling where demixing or phase separation is 
one of many unknowns. Thus, laboratory experiments 
are needed to investigate these extreme conditions. This 
is even more true for such complex, but very interesting 
processes as the formation of pure carbon from methane 
under the high pressures found in Neptune [8|. 

Driven by the progress in inertial confinement fusion 
[9|, [lOJ, mixing receives currently also much interest in 
laboratory settings. As the boundary between the fuel 
and the ablator is hydrodynamically unstable during the 
compression phase, fuel and shell material (beryllium or 
carbon) mix. The degree of mixing strongly affects the 
performance of the target and can even prevent sufficient 



burn or ignition at all Uj, ll2|. However, experimental 



techniques to probe mixing under the extreme conditions 
in fusion targets prove to be very difficult and only a few 
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spectroscopic investigations exist 

however that present models for the mixing during the 

compression underestimate the experimental results. 

In this Letter, we develop the theoretical framework 
for the application of elastic x-ray scattering as a tool 



to investigate mixtures and the high-pressure miscibihty 
of materials. X-ray Thomson scattering has been shown 
to robustly deliver basic plasma parameters like density, 
temperature, and ion charge state as well a s dy namic and 
structural properties of simple materials 15l4l7|. First 
experiments with two-component systems like plastics 
(CH) or lithium hydrate (LiH) have been performed as 
well [181421 1 . Hitherto the theoretical description of the 
spectrum was however based on an ion structure that was 
generated from a one-component calculation via |18| 
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Here, Sah denote the partial structure factors, n^ and Za 
are the density and the charge of the ions of species a. 
The structure for the effective one-component system, i.e. 
S\^ , is calculated using an average ion charge state de- 
fined by Z — ^^ naZal X^a "■a- Such a treatment is only 
exact in the limit of weakly coupled plasmas that can 
be described by the random phase approximation (RPA) 
[22l |. It is however unable to describe the highly nonlin- 
ear effects in the structure of strongly coupled ions [23| . 
Therefore, to apply x-ray scattering as a reliable diag- 
nostic method for dense matter requires a theoretical de- 
scription that considers the full microscopic structure in 
the materials tested including the nonlinear interplay be- 
tween different highly correlated ion species. In contrast 
to the approach of Ref. [l8(, the new approach should 
also allow for Z-dependent screening clouds for different 
ion species. 

Here, we derive a full multi-component description of 
the x-ray scattering signal. The partial structure factors 
required are obtained via density functional molecular 
dynamics simulations (details in Refs. [2J,|25|). We show 
that the elastic scattering feature is very sensitive to the 
ratio of the different elements in the scattering volume. 
Thus, it can be used as a probe for the degree of mixing 
in strongly compressed samples. Moreover, we predict 
considerable differences in the scattering strength from 
microscopically mixed systems and matter consisting of 
two phases. The differences are particularly pronounced 



for hydrogen-helium mixtures under conditions found in 
the interior of Jupiter. Thus, elastic x-ray scattering can 
be used to investigate demixing of hydrogen and helium 
under planetary conditions. 

We briefly sketch our theoretical description of x-ray 
Thomson scattering that is developed along the ideas of 
Chihara [26|, |27|, but fully generalized to plasmas with 
multiple ion species. The influence of the ions is however 
quite indirect as x-rays are scattered by electron density 
fluctuations. The scattered intensity is thus proportional 
to the total electron structure factor S'g°*(k, w), where a; 
and k are the frequency and wave number shifts of the 
photon, respectively [28|. The structure factor can be 
expressed as the Fourier-transform of the intermediate 
scattering function with respect to time. The latter is 
the correlation function of electron densities 



Fr(k,t) = (gr(k,i)er(-k,o)). 



(2) 



Now we split the total electron density into free elec- 
trons and contributions of core electrons associated with 
N different ion species: gl°\k, t) = Ef qW^^ t)+QHK t). 
Applying this decomposition in the quadratic form ([2]), 
we obtain N'^ bound-bound terms, N bound-free terms 
and one free-free term. All of these terms can be treated 
as in the case of just one ion species. Special care must 
be however given to core electrons belonging to different 
ions species X]a^f)(^a(k, ^^^("k, 0)) as it defines a dis- 
tinct ion-ion scattering function without core excitations. 
Evaluating the different density correlations, we obtain 
for the total electron structure factor (in the following, 
the fc-dependence is dropped for simplicity) 



a 

a -^ 

+ y^ ^/XgXb fafb Sabjio) ■ (3) 

a, 6 

Here, we introduced the concentrations Xa = Ua/ X^a "■a- 
Z'^ is the number of core electrons bound to ions of the 
species a. The first term describes correlations between 
two free electrons. The next term accounts for free-bound 
correlations where fa is the atomic/ionic form factor of 
bound states of component a. The second line contains 
self-contributions, i.e., internal excitations and bound- 
free transitions. Except the summation over species, it is 
unchanged from its usual form. The last term describes 
correlations between two bound electrons. 

The main problem using the expression ^ is the fact 
that all partial structure factors are connected. This 
becomes particularly clear when considering See which 
contains correlations between two screening clouds and, 
thus, also ionic properties. Indeed, the structure factors 
used in Eq. © form a set of \N{N + l) equations poj . 



This set can be conveniently written in matrix form and 
then inverted. Such a procedure yields 

ZSeaij^) = Xa Qa Saa{i^) +^Xb Qb Sab{i^) (4) 

for the free-bound structure factor. Here, Sab denotes the 
partial ion-ion structure factors. The correlations of the 
free electrons to the ions are contained in the screening 
fmiction qa{k) which are defined via 



<la{k) 



neCea{k)xl{k) 
l-neCee(fc)x°(fc)' 



(5) 



where Cee and Cea are the direct electron-electron and 
electron-ion correlation functions, respectively. In lowest 
order, these are given by the respective potentials. Y°(fc) 
denotes the density response of a free electron gas [29| . 
For the free-free structure factor, one obtains 



ZSee{^) 



a.b 



y/XaXb qaQb Sabil^) + Seei^) , (6) 



which now separates electron-ion correlations from the 
structure factor of the free electron gas S^^{k,uj) that 
characterizes the kinetically free electrons in the system. 
The results can be summarized in the total electron 
structure factor of the form 

SlTioj) = J2 V^^ [fa + qa][fb + Qb] Sabiio) + Z S°e{uj) 
a,b 

+ J2 K Xa f dw' ~ST{oj - uj')S!{u;') . (7) 

The first term describes quasi-elastic scattering at bound 
electrons and the screening clouds associated to different 
ion species. In this contribution, the full ionic structure, 
expressed by the partial structure factors Sabi influences 
the scattering spectrum. As the ion motion cannot be 
resolved in current laser experiments, the ion structure 
can be treated statically: Sabik,io) — Sab{k)x6{uj). The 
static structure factors Sabik) can be obtained by means 
of classical hypernetted chain (HNC) calculations [23[ or 
numerical simulations |25l . |30[ • The second term contains 
the full dynamic response of the free electron gas which 
can be in lowest order described by the random phase 
approximation. Extensions include electron-ion collisions 
and local field corrections 3l|. The last term describes 
contributions due to the the excitation or ionization of 
bound electrons by x-rays. 

The inelastic contributions due to scattering at free 
electrons, bound-free transitions, and internal excitations 
are unchanged from a description for systems with one 
ion species [18[. We thus concentrate here on elastic 



scattering, that is the first term of Eq. ([7]), which high- 
lights the mutual correlations between the different ion 
species. To calculate the weight of the Rayleigh peak, i.e.. 
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FIG. 1: Comparison of different treatments for tlie weigiit of 
tlie Rayleigli pealc Wij(fc) (see text). Considered is a CH- 
plasma with nc ~ nn ~ 2.3 x 10^^ cm~'^ and a temperature of 
T = leV. The ion charge states are Zc~4: and Zh = 1. Both 
IC and partial structure factors are calculated via HNC [2a |. 



WRik)=j:^,^II^[faik)+qam[fb{k)+qbik)]Sabik), 

we use structure factors from ab initio simulations (DFT- 
MD) (2 a or solutions of the hypernetted chain equations 
(HNC) [23|. The screening functions are used in linear 
response to a Coulomb field, i.e., qa{k) — ZaK?/{n'^ + k'^)^ 
where k is the inverse electron screening length. The 



form factors fa are taken from isolated atoms/ions [32 1. 
First we compare the full multi-component formula ([7]) 
and the approximate treatment of Ref. [18] which is based 
on the structure of an average ion component via Eq. ([T]) . 
Both approaches agree for weakly coupled systems, but 
large differences occur for strongly coupled systems as 
shown in Fig. [T] {Tu « 80). For both the screened and un- 
screened ion systems, the maximum of the Rayleigh peak 
is shifted. This reflects the fact that the higher charged 
carbon ions imprint their structure into the proton sub- 
system which cannot be described by the one-component 
structure. Moreover, the mutual screening of the ions is 
neglected in the reduced model which results in a strongly 
underestimated Rayleigh peak at small k. The screened 
interactions are more realistic 171, but the results for un- 



screened ions demonstrate that the reduced model may 
even predict negative Wr (especially, for small k). For 
strongly coupled multi-component systems, the analysis 
should thus be based on the new expression ([7]). 

Let us now turn to the application of the new multi- 
component description for the mixing of beryllium and 
hydrogen as it occurs during inertial fusion experiments. 
Due to the strong drive, the initially well-defined inter- 
face between the two materials experiences a Rayleigh- 
Taylor instability. At this stage, a volume element close 
to the original boundary will contain two fluids. However, 
these fluids consist of either pure beryllium or hydrogen 
and the system is made of two distinct phases. Later 
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FIG. 2: (color online) The upper panel shows the weight of 
the Rayleigh peak, WR{k), for real microscopic mixtures of 
beryllium and hydrogen compared to systems that contain 
both materials in two pure phases. The concentrations of 
beryllium are (from top to bottom): 100%, 80%, 50%, 20%, 
and 0%. The temperature of the system is T = 13 eV. The 
total densities were arranged to be in pressure equilibrium 
with pure beryllium at riBe = 3.707 x lO^'^cm"'^. The lower 
panel displays the partial beryllium-beryllium structure factor 
for the same concentrations. 



both materials will microscopically mix due to diffusion. 

Figure [2] demonstrates that the degree of mixing in 
dense beryllium-hydrogen systems can be investigated by 
measuring the strength of elastic x-ray scattering. To 
allow for experimental verification, we have taken the 
density and the temperature from a recent experiment 
on shocked pure beryllium 33[. Behind the shock front, 
all possible mixtures will be compressed until they are 
in pressure equilibrium with the driven beryllium. Thus, 
we have adjusted the total densities of the mixtures until 
the pressure matches the one of the pure beryllium. This 
procedure yields, e.g., a density of tih = 8x 10^'^ cm^^^ for 
the case of 100% hydrogen. 

As beryllium scatters more efficiently than hydrogen, 
the elastic feature is a strong function of the mixing ratio. 
Interestingly, both microscopic mixtures and two-phase 
systems scatter very similar at small k. Here, the signal 
is mainly determined by the partial beryllium-beryllium 
structure factor and we can use these data to determine 
the concentrations. However, significant differences arise 
at the maximum of the elastic scattering peak around 
A: = 5 A~^ that can be used to distinguish between hydro- 
dynamic and diffusive mixing. 
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FIG. 3: (color online) Strength of the elastic scattering peak 
for pure hydrogen, pure helium, and a mixture with xyl = 0.924 
and a::Hc = 0.076 at a temperature of r = 5000K (upper panel). 
All systems have an electron density of rie =4.7x 10 cm~ . 
The lower panel shows the partial structure factors of the 
three systems above. 



In the next example, we investigate a hydrogen-helium 
mixture under conditions as found in the interior of giant 
gas planets [xhc = "-Ho/('t-Hc + "-h) = 0.0756]. Density 
and temperature yield a system in the atomic/molecular 
phase (no ionization). Fig. [3] shows that the strength 
of the elastic scattering feature displays large differences 
between the pure substances and the mixture. Indeed, 
we find very distinct scattering features for the different 
systems: pure helium scatters by far the most effective 
and displays a peak around fc = 4 A^^, the mixture shows 
a monotonically decreasing shape, and Wjj(fc) for pure 
hydrogen is almost featureless. Due to the mixing ratio, 
the two phase system is dominated by the properties of 
hydrogen. 

The reason for this behavior lies not only in the fact 
that helium atoms scatter the x-rays more efficiently. The 
different structure factors play also an important role. 
These structure factors were calculated from ah initio 
quantum simulations (see Ref. [2J]) which allows us to 
treat neutral systems as well. Pure helium has a strong 
peak in the structure factor which is also clearly visible in 
the elastic scattering peak. This feature disappears in the 
mixture as the helium concentration is low. However, the 



small fraction of helium mitigates the correlations within 
the hydrogen subsystem. Thus, small hydrogen-helium 
correlations lead to a shape of Wr that is dominated by 
the decreasing atomic form factor. Similar results are 
also found for denser and hotter systems where hydrogen 
and helium are ionized (not shown). 

In conclusion, we have derived a novel description for 
x-ray scattering in systems with multiple atomic or ion 
species that take all the mutual correlations into account. 
The results show that x-ray scattering is a powerful tool 
to investigate dense mixtures and the mixing/demixing 
of materials under extreme conditions as they are found 
during inertial confinement fusion and in the interior of 
of old stars and giant planets. 
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